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We examine the level of enhancement that can be achieved in the ZZ and 77 channels for a 
two-Higgs-doublet model Higgs boson (either the light h or the heavy H) with mass near 125 GeV 
after imposing all constraints from LEP data, B physics, precision electroweak data, vacuum sta- 
bility, unitarity and perturbativity. The latter constraints restrict substantially the possibilities for 
enhancing the gg — > h — > 77 or gg — > H — > 77 signal relative to that for the SM Higgs, /ism- 
Further, we find that a significant enhancement of the gg — > h — > 77 or gg — s> H — s> 77 signal in 
Type II models is possible only if the gg — > h — > ZZ or gg — > H — > ZZ mode is even more enhanced, 
a situation disfavored by current data. In contrast, in the Type I model one can achieve enhanced 
rates in the 77 final state for the h while having the ZZ mode at or below the SM rate — the largest 
[gg — > h — > 77]/[g<? — > /ism77] ratio found is of order ~ 1.3 when the two Higgs doublet vacuum 
expectation ratio is tan j3 — 4 or 20 and the charged Higgs boson has its minimal LEP-allowed value 
of m H ± = 90 GeV. 
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Data from the ATLAS and CMS collaborations [TJ |5] provide an essentially 5a signal for a Higgs-like resonance with 
mass of order 123-128 GeV. In the 77 final state, the ATLAS and CMS gluon fusion induced rates are significantly 
enhanced relative to the the Standard Model (SM) prediction. In the gg — >• h — >• ZZ —> At channel, the ATLAS and 
CMS signals are not far from the SM expectation, with ATLAS showing a slight enhancement and CMS showing 
a slight suppression. In the bb, t + t~ and WW — > tvlv channels, the central value ATLAS rates are somewhat 
suppressed relative to the SM prediction but error bars are very large. The CMS signals in these latter channels 
are also somewhat suppressed and lie at least la below the SM prediction — no signal being observed in the t + t~ 
channel. Meanwhile, the CDF and DO experiments have announced new results [3] that support the ~ 125 GeV Higgs 
signal and suggest an enhancement relative to the SM of the VF+Higgs with Higgs— > bb rate by a factor of 2 ± 0.6. 

Enhancements with respect to the SM in the 77 channel are generically possible in two-Higgs-doublet models 
(2HDM) of Type-I and Type-II as explored in |4j[5]. However, these papers do not make clear what level of enhance- 
ment is possible after all constraints from B physics and LEP data (B/LEP), precision electroweak data, unitarity 
and perturbativity are imposed. In this paper, we impose all such constraints and determine the maximum possible 
enhancement. We employ a full 1-loop amplitude for Higgs— > 77 without neglecting any contributions from possible 
states in the loop. We examine correlations with other channels. We also consider cases of degenerate scalar masses 
at ~ 125 GeV 0. 
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I. INTRODUCTION 
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II. 2HDM MODELS 



The general Higgs sector potential employed is 

V =m 2 1 $|$i + m! 2 3>2^2 - m 2 2 <f>|$ 2 + h.c. 

+ iAi ($I$i) 2 + ^A 2 ($|$ 2 ) 2 + A 3 (*t $1 ^ ($t $2 ) + A4 ^t$ 2 ) ^t $i ^j 



(1) 



^A 5 ($1$ 2 ) 2 + [a 6 + A 7 ($|* 2 )] ($1*2 ) li.c 



where, to avoid explicit CP violation in the Higgs sector, all Aj and m\ 2 are assumed to be real. We choose a basis 
in which 



lib \ - — ( 
[ 1} ~ V2 \cosf3 



(*2> = — ( r° 



where u = (v^Gf) -1 ^ 2 ~ 246 GcV. By convention < /3 < n/2 is chosen. For real parameters, the phase £ could 
still be non-zero if the vacuum breaks CP spontaneously. We avoid parameter choices for which this happens and 
take £ = 0. Then, we define 



with vi = vcos/3 and t> 2 
pseudoscalar state is 



The physical scalars are: 



v sin /3. 



The neutral Goldstone boson is G° 
A = — ?7i sin /3 + r] 2 cos /3 . 



1,2 (2) 
771 cos /3 + rj2 sin /3 while the physical 



h — —pi sin a + P2 cos a, H — pi cos a + p 2 sin a . 



(3) 



(4) 



Without loss of generality, one can assume that the mixing angle a varies between — 7r/2 and w/2. We choose our 
independent variables to be tan/3 and sin a, which are single valued in the allowed ranges. 

We adopt the code 2HDMC [7] for numerical calculations. A number of different input sets can be used in the 
2HDMC context. We have chosen to use the "physical basis" in which the inputs are the physical Higgs masses 
rnft, m^4, mjy±), the vacuum expectation value ratio (tan/3), and the CP-even Higgs mixing angle, a, supple- 
mented by m\ 2 - The additional parameters A6 and A7 are assumed to be zero as a result of a Z 2 symmetry being 
imposed on the dim 4 operators under which Hi — > Hi and Hi — > — i? 2 . m\ 2 7^ is still allowed as a "soft" breaking 
of the Z2 symmetry. With the above inputs, Aj.,2,3,4,5 as well as m 2 ^ and m| 2 are determined (the latter two via the 
minimization conditions for a minimum of the vacuum) [5] . 

In this paper we discuss the Type I and Type II 2HDM models, that are defined by the fermion coupling patterns 
as specified in Table [I] — for more details see [5] . 





Type I 


Type II 


Higgs 


up quarks 


down quarks 


leptons 


up quarks 


down quarks 


leptons 


h 


cos a/ sin /3 


cos a/ sin j3 


cos a/ sin /3 


cos a/ sin /3 


—sin a/ cos j3 


—sin a/ cos j3 


H 


sin a/ sin (3 


sin a/ sin /3 


sin a/ sin /3 


sin a/ sin /3 


cos a/ cos /3 


cos a/ cos /3 


A 


cot /3 


-cot/3 


-cot/3 


cot/3 


tan /3 


tan /3 



TABLE I: Fermionic couplings C ^ normalized to their SM values for the Type I and Type II two-Higgs-doublet models. 



III. SETUP OF THE ANALYSIS 



The 2HDMC code implements precision electroweak constraints (denoted STU) and limits coming from requiring 
vacuum stability, unitarity and coupling-constant perturbativity (denoted jointly as SUP). The SUP constraints are 
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particularly crucial in limiting the level of enhancement of the gg —¥ h — > 77 channel, which is our main focus. For all 
our scans, we have supplemented the 2HDMC code by including the B/LEP constraints. For the LEP data we adopt 
upper limits on er(e + e~ — > Zh/H) and a(e + e~ — > Ah/H) from [TU] and [TT], respectively. 1 Regarding B physics, 
the constraints imposed are those from BR(£> S — > X s j), Rj,, AMb,, £k, BR(£? + — > t + v t ) and BR(£> + — > Dt + v t ). 
The most important implications of these results are to place a lower bound on m H ± as a function of tan/3 as shown 
in Fig. 15 of [12] in the case of the Type II model and to place a lower bound on tan j3 as a function of m H ± as shown 
in Fig. 18 of [H]. 

While looking for an enhancement of the signal in the 77 channel we also computed the extra Higgs-sector contri- 
butions to the anomalous magnetic moment of the muon, = (g^ — 2)/2. Since the experimentally measured value, 
dju = (1165920.80 ± 0.63) x 10 -9 [13], differs by ~ 3cr from its SM value it is important to check correlations between 
Sa^ = a M — a^ M and the signal in the 77 channel. Given the B/LEP, STU and SUP constraints, it turns out that 
one-loop contributions within the 2HDM are small and negligible, and the leading contribution is that known as the 
Barr-Zee diagram |14] which emerges at the two-loop level. For completeness we include also sub-leading contributions, 
see [7]. Since the overall ~ 3er discrepancy between the experimental and theoretical SM values could still be due to 
fluctuations (the world average is based mainly on the E821 result [TS] with uncertainties dominated by statistics) or 
underestimates of the theoretical uncertainties, we do not use the a M measurement as an experimental constraint on 
the models we discuss. However, in tables presented hereafter we do show (in the very last column in units of 10~ n ) 
8a^, the judgment as to whether <5a M is acceptable being left to the reader. In fact, for all parameter choices yielding 
an enhanced Higgs to two-photon rate the extra contributions to are very small and the a M discrepancy is not 
resolved. 

Coupling constant perturbativity is defined as in 2HDMC by the requirement that all self-couplings among the 
Higgs-boson mass eigenstates be smaller than 4-7T. For the scenarios we consider, this becomes an important constraint 
on Ai. By the vacuum stability constraint we mean positivity of the potential in all directions at asymptotically large 
field strength. Unitarity is the requirement that the multi-channel Higgs scattering matrix have a largest eigenvalue 
below the unitarity limit using the analysis of [16) . 

For an individual Higgs, denoted hi (where hi — h,H,A are the choices) we compute the ratio of the gg or WW- 
fusion (VBF) induced Higgs cross section times the Higgs branching ratio to a given final state, X, relative to the 
corresponding value for the SM Higgs boson as follows: 

R h i(Y\ = (r h i\ 2 B R( h i^ x ) nhi (v\ = (r< hi \2 BR(hj^X) 

where hg M is the SM Higgs boson with m hsM = m h . and Cg* , C^ w are the ratios of the gg — ¥ hi, WW — > hi couplings 
{Pww being zero at tree level) to those for the SM, respectively. Note that the corresponding ratio for V* — > Vhi 
{V = W, Z) with hi -> X is equal to R^ BF (X), given that kinematic factors cancel out of all these ratios and that 
these ratios are computed in a self-consistent manner (that is, treating radiative corrections for the SM Higgs boson 
in the same manner as for the 2HDM Higgs bosons). When considering cases where more than one hi has mass of 
125 GeV [§], we sum the different R hi for the production/decay channel of interest. 
We have performed five scans over the parameter space with the range of variation specified in Table [H] 





scenario I 


scenario II 


scenario III 


scenario IV 


scenario V 


m h [GeV] 


125 


{10,..., 124.9} 


125 


125 


{10,..., 124.9} 


m H [GeV] 


125+{0.1,. . . ,1000} 


125 


125.1 


125+{0.1,...,1000} 


125 


m A [GeV] 


{10,..., 1000} 


{10,..., 1000} 


{10,..., 1000} 


125.1 


125.1 


m H ± [GeV] 


1500 (tan^=0.5); 800 (tan ,0=1); 250,350 (tan/3=2); 90,150,250,350 (tan/3 > 2) for Type I 
600 (tan/3=0.5); 500 (tan^=l); 340 (tan/3 = 2); 320 (tan/3 > 2) for Type II 


tan p 


{0.5, ...,20} 


sin a 


{-1, ...,1} 


m'i 2 [GeV] 


{-1000 2 ,. ..,1000 2 } 



TABLE II: Range of parameters adopted in the scans. The values of m H ± are bounded from below by the constraints from B 
physics, see Fig. 15 and Fig. 18 of [12] for the Type II and Type I models, respectively. 



We have modified the subroutine in 2HDMC that calculates the Higgs boson decays to 77 and also the part of the code relevant for 
QCD corrections to the qq final state. 
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A. The rrih = 125 GeV or tuh = 125 GeV scenarios: 



Let us begin by discussing the case in which the h has mass = 125 GeV while scanning over the masses of the 
other Higgs eigenstates (degenerate cases are discussed below) . In Fig. [I] we plot the maximum value achieved for the 
ratio Rg g (77) as a function of tan f3 after scanning over all other input parameters (as specified earlier) , in particular 
sin a. These maximum values are plotted prior to imposing any constraints and after imposing various combinations 
of the constraints outlined earlier with point notation as specified in the figure legend. We observe that for most 
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FIG. 1: The top two plots show the maximum Rggi^fy) values in the Type I (left) and Type II (right) models for mu = 125 GeV 
as a function of tan/3 after imposing various constraints — see figure legend. Corresponding Rg g (ZZ) and Rg g (bb) are shown 
in the middle and lower panels. Disappearance of a point after imposing a given constraint set means that the point did not 
satisfy that set of constraints. In the case of boxes and circles, if a given point satisfies subsequent constraints then the resulting 
color is chosen according to the color ordering shown in the legend, the same pattern is adopted in the remaining plots. 
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FIG. 2: Rggi'Y'y) is plotted for mj, = 125 GeV as a function of sina for a sequence of tan/3 values. Different constraint 
combinations are considered and the different curves of given type correspond to a variety of other input parameters. The 
upper plots are for the Type I model and the lower plots are for the Type II model. Different colors indicate different tan/3 
values. 



values of tan/3 the B/LEP and STU precision electroweak constraints, both individually and in combination, leave 
the maximum Rggill) unchanged relative to a full scan over all of parameter space. In contrast, the SUP constraints 
greatly reduce the maximum value of Rg g (77) that can be achieved and that value is left unchanged when B /LEP and 
STU constraints are imposed in addition. Remarkably, in the Type I model maximum R gg (^^) values much above 1.3 
are not possible, with values close to 1 being more typical for most tan/3 values, whereas maximum Rggij'f) values 
in the range of 2 — 3 are possible for 2 < tan/3 < 7 and tan/3 = 20 in the Type II model. 

One can get a feeling for how the different constraints impact Rggij'f) by plotting this quantity as a function of 
sina at fixed tan/3 for different constraint combinations and a selection of different other input parameters. As shown 
in Fig. 2 B^J'yy) typically has a maximum as sina is varied but the height of this maximum depends very much on 
the const rain ts i mpo sed as there is also variation with the other input parameters. 



Tables 



III 



and 



IV 



display the full set of input parameters corresponding to the maximal i?g S (77) values at each 
tan /3 for models oFType I and Type II, respectively. It is important to notice that in the Type II model, the value of 
Rgg(ZZ) corresponding to the parameters that maximize Rg g (-fj) is typically large, ~ 3. In fact, as discussed shortly, 
Rg JZZ) > Rggill) whenever ii^ (77) is even modestly enhanced. Since, experimentally, the observed signal for ZZ 
is substantially less enhanced than that for 77, the Type II models seem to be disfavored. We note that, in contrast, 
Rgg(ll)/Rgg( zz ) Z 1 when Rg g hl) > 1 is fairly typical of the MSSM model (which has a Type II Higgs sector), 
especially with full or partial GUT scale unification for the soft-SUSY-breaking parameters, see for example [17] . In 
such scenarios the primary modification to the 77 rate relative to the SM is due to the light stop loop contribution 
to the /177 coupling (which enters with the same sign as the W loop and has a color factor enhancement) which 
enhances BR(/i — > 77). Note that the stop loop contribution to the hgg production coupling is the same for both the 
ZZ and 77 final states. In the absence of GUT scale unification, there are many other potentially significant loops 
contributing to an increase in the /177 coupling, the most important being the light chargino loop and the light stau 
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0.98 
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400 
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0.99 
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0.99 
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90 
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-1.7 


5.0 


0.98 


0.98 


1.06 


0.90 


0.91 


0.98 
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400 


150 


100 


-0.0 


0.01 


-1.6 


7.0 


1.04 


0.99 


0.98 


1.06 


1.01 


0.99 


135 


500 


90 


50 


-0.2 


0.02 


-1.6 


10.0 


0.90 


0.81 


0.74 


0.99 


0.89 


0.81 


175 


500 


150 


50 


-0.5 


0.04 


-1.5 


15.0 


0.46 


0.59 


0.66 


0.41 


0.53 


0.59 


225 


400 


350 


50 


0.6 


-0.11 


-1.4 


20.0 


1.31 


1.00 


1.00 


1.30 


0.99 


1.00 
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-0.0 


0.13 


-1.5 



TABLE III: Table of maximum Rggijj) values for the Type I 2HDM with rrih = 125 GeV and associated R values for other 
initial and/or final states. The input parameters that give the maximal R^ g {'yy) value are also tabulated. 
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10.0 


0.34 
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0.28 
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0.87 


225 


200 


320 


50 
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-0.04 


2.4 



TABLE IV: Table of maximum -Rg ff (77) values for the Type II 2HDM with mt = 125 GeV and associated R values for other 
initial and/or final states. The input parameters that give the maximal Rggl'yy) value are also tabulated. 



loop, as studied for example in [T5] . 

Corresponding results for the H are presented for the Type I and Type II models in Tables [V] and |VI| respectively. 
In the case of the Type I model, an enhanced gluon fusion rate in the 77 final state does not seem to be possible after 
imposing the SUP constraints, whereas maximal enhancements of order ^^(77) ~ 2.8 are quite typical for the Type 
II model, albeit with even larger R^ g (ZZ). Again, in the case of the Type II model #^(77) / 'R^ g (Z Z) < 1 applies 
more generally whenever iij^ (77) is significantly enhanced. 

That an enhanced 77 rate, e.g. R gg H (77) > 1.2, leads to R gg H (77) / R gg H (Z Z) < 1 in Type II models is illustrated 
by the plots of Fig.pl We again emphasize that this is to be contrasted with the Type I model for which Rgg^j) > 1.2 
implies R gg (^^) / R]^ 1 (Z Z) > 1, see Fig.jiJ in better agreement with current data. (For the Type I model, R^ g (-f-f) > 1 
is not possible after imposing the SUP constraints.) 

It is interesting to understand the mechanism behind the enhancement of R gg H (ZZ) that seems to be an inevitable 

result within the Type II model if R gg H (jj) is large. Let us define r as the ratio of 77 over ZZ production rates for 
a scalar s (either h or H). Then it is easy to see that 

^g g (77) = T(.s -» 77)/F(fr sm -> 77) (6) 
Ts R s gg {ZZ) T(s -> ZZ)/T{h sm -+ZZ)' 1 ' 

For the decay mode s — > ZZ* , the tree level amplitude is present and dominant so that the numerator simply reduces 
to [C s zz ) 2 . For the decay mode s — > 77, there is no tree level contribution — the 577 coupling first arises at the 
one-loop level with the t-loop, WToop and _ff ± -loop being the important contributions. As a result, the denominator 
can be written as 

r(* -> 77) = ( CtvwA^ 1 - C s ti A s ; M + A H ± 
r(/w l ^77) I A s w M -Af M 

where and C^ vw are the sti and sWW couplings normalized to those of the hsu., while A^ 1 an d Af M are the 
lU-loop and t-loop amplitudes, respectively, for the /ism- Finally, An± is the iJ ± -loop amplitude in the 2HDM; since 
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FIG. 3: For the h (top) and (bottom) we plot the gluon fusion induced "/"//ZZ ratio as a function of Rggi'J'y) for the Type 
II 2HDM. 



it is very small in the Type II model, it can be neglected. Thus, 

{Cww) \ Aw c ww 1 I w C'ww 1 



(T s \Z ASM _ ASM ASM _ ASM 



(8) 



where C s zz = C ww in any doublets+singlets models. Note that when the i-loop contribution is negligible then 
r s — > 1. It is easy to see that r s < 1 if the following inequality is satisfied 



ns asm 

'-'WW 



When C^/C ww is outside of the above interval then r s > 1. If s is the lighter scalar /i then C^/C^ w 
cosa/[sin/3sin(/3 — a)] implying < 1 when 



^sin/Jsin^-a) <2 ^p _1 - 9 ' (10) 
while for s — H, C^/C ww — sina/[sin/3cos(/3 — a)] and we obtain r# < I for 

1< ^ <2 ^--l~9 (11) 

In the case s = h, Rggij'f) is maximized by suppressing the ft. total width, which corresponds to chosing a so as to 
minimize the hbb coupling, i.e. a ~ 0, resulting in Ci/(7^r W ~ 1/ sin 2 /? > I (and < 5 for tan/3 > 0.5). Consequently 



Th < 1, Sts observed in Table IV The argument is similar in the case of the H: this time the Hbb coupling is chosen to 
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FIG. 4: For the /i (top) and H (bottom) we plot the gluon fusion induced 'yy/ZZ ratio as a function of Rggi'J'y) for the Type 
I 2HDM1. 



be small (equivalent toa~ ±7r/2) in order to minimize the H total width and therefore maximize i?^(77), with the 
result that once again /C^w ~ 1/ sin 2 (3 > 1, yielding rjj < 1. These analytic results explain why large R s gg {^) 
is correlated with even larger R gg (ZZ) in Type II 2HDMs. 

In Fig. [5] we plot contours of r s and of R s gg {ll) m the Type II model. It is seen from the left panel that if tan (3 
is large then only small a's will maximize ^^(77). And, in that region, is always less than 1. Note that the 
Rggill) > 1 region shrinks for large tan/3, so the the values of a preferred for large R gg (ll) converge to when 
/3 —> it/2. For the case of s = H, the right panel shows that when tan/3 is large then only vertical bands of a 
corresponding to values close to ±?t/2 are allowed if R gg (jj) > 1. From the plots, we see that g ( 77) > 1 could be 
consistent with r s > 1 only if tan/3 < 1, which explains the pattern observed in Tables IV and VI and FigjSj Note, 
however, that small tan /3 is disfavored by B-physics as it enhances the H + tb coupling too much, see for example |12) . 

Once again, we emphasize that a substantial enhancement of the 77 rate is possible for the h in Type I models 
without enhancing the ZZ rate. In particular, from Table III we see that the enhancement in the 77 channel is ~ 1.3 
(for both gg fusion and VBF) for tan (3 — 4 and 20 while other final states, in particular ZZ, have close to SM rates. 
The table also shows that this maximum is achieved for sin a ~ 0. Thus, (3 ~ 7r/2 and cos a ~ 1 yielding SM-like 
coupling of the h to quarks (see Table |l| and vector bosons. It turns out that in these cases the total enhancement, 
~ 30%, is provided by the charged Higgs boson loop contribution to the 77-coupling. In these same cases, the mass 
of the heavier Higgs boson is run — 225 GeV. As such a mass is within the reach of the LHC, it is important to make 
sure that the H cannot be detected (at least with the current data set). It is easy to see that indeed this is the case. 
Since guzz oc cos(/3 — a) and g Hb i Ht j oc sin a one finds that the H decouples from both vector bosons and fcrmions 
given that a ~ and j3 ~ tt/2. The A will also be difficult to detect since it has no tree-level WW, ZZ coupling and the 



III 



we 



Abb, Att couplings, being proportional to cot /3, will be quite suppressed, especially at tan /3 = 20. From Table 
observe that for tan/3 = 4 and 20 the corresponding charged Higgs is light, m H ± = 90 GeV, i.e. as small as allowed by 
LEP2 direct searches in e + e~ — > H + H~ . Searches for a light are underway at the LHC along the lines described 
in [TU]. The most promising H ± production and decay process is pp —> tt — >• H ± bW T b —> rvbbq'q. According to Fig. 3 
of [19] , for the Type I model, the region of tan /3 < 6 — 7 for m#± ~ 90 GeV could be efficiently explored at the 14 TeV 
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TABLE V: Table of maximum ^^(77) values for the Type I 2HDM with mn = 125 GeV and associated R values for other 
initial and/or final states. The input parameters that give the maximal Rfgi'yy) value are also tabulated. 
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TABLE VI: Table of maximum Rfgijj) values for the Type II 2HDM with ran = 125 GeV and associated R values for other 
initial and/or final states. The input parameters that give the maximal Rfgi'y'y) value are also tabulated. 




FIG. 5: In the left panel we show contour plots (with shadowing) in the (f3, a) space for Th with superimposed red contours 
for -Rg 9 (77). The right panel show similarly ru with ^^(77). Red numbers give constant values of _Rg g (77) (Rfg("f"f)) while 
black ones show constant values of ru (th)- The white region correspond to r s > 10.75. 



LHC by ATLAS even at the integrated luminosity of 10 fb 1 - for more details see [19]. The existing LHC bounds on 
BR(t H+b) obtained assuming BR(i? ± t±v t ) = 1 are only moderately restrictive: 5% - 1% [ID] (4% - 2%)[H] 
for masses of the charged Higgs boson m H ± = 90(80) — 160 GeV in the case of ATLAS (CMS), respectively. These 
bounds are weakened in the Type I model where Si?(i7 ± — > r v T ) ~ 0.7. Since BR(t — > H + b) ~ 1/ tan 2 f3, large 
tan/3 suppresses BR(t — > H + b). Indeed, it is easy to verify that for m H ± = 90 GeV BR{t — > H + b) is ~ 3.8% and 
~ 0.15% for tan/3 = 4 and tan/3 = 20, respectively. So, a charged Higgs yielding enhanced h —> 77 rates in gg fusion 
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and VBF is still completely consistent with current data. 



B. rriA ~ 125 GeV and rrih, = 125 GeV or rriH = 125 GeV scenario 



The signal at 125 GeV cannot be pure A since the A does not couple to ZZ, a final state that is definitely present 
at 125 GeV. However, one can imagine that the CP-even h or H and the A both have mass of 125 GeV and that the 
net 77 rate gets substantial contributions from both the h or H and the A while only the former contributes to the 
ZZ rate. These possibilities are explored in Figs. [6] and [7j from which we observe that an enhanced 77 rate is only 
possible for the rrih — rriA — 125 GeV choice . De tails for this case appear in Table VII 
For the Type I model, we see from Table 



VII 



that i?g ff (77) is significantly enhanced only for the same tan/3 = 4 
and tan j3 = 20 values as in the case of having (only) rrih = 125 GeV and that the pseudoscalar contribution i?^ g (77) 
turns out to be tiny. However, the contribution to the bb final state from the A can be substantial. Given that 
the top loop dominates both the Agg and hgg coupling one finds (C A g /Cg g ) 2 ~ (3/2) 2 (cos/3/ cos a) 2 , where we used 
C^/C^ — cos/3/cosa from Table |l] and the m,h,A 2mt fermionic loop ratio of A/h = 3/2. As a result, the A can 

contribute even more to the bb final state rate than the h if tan/3 is small. This (unw anted) contribution to the bb 

for tan /3 — 2 — 4. In the end, 



VII 



final state from A production is apparent from the results for R^ A {bb) in Table 
only tan/3 = 20 yields both an enhanced 77 rate, Rg^ ax {jj) — 1.31, and SM-like rates for the ZZ and bb final 
states, R^+ A (ZZ) = R h + A {bb) = 1. For this case /3 ~ it/2 and a = implying that the h couples to fermions and 
gauge bosons like a SM Higgs boson and the enhancement of -Rg S h ,nax(77) * s due exclusively to the charged Higgs loop 
contribution to the 77 couplings. 

For the Type II model, see Table VIII the pseudoscalar contribution Rfgi'f'f) is also (as for the Type I model) 
negligible. Thus, the enhancement of -Rg^' 4 (77) is essentially the same as that for i?Jj g (77) for the case when only 
nih = 125 GeV, reaching maximum values of order 2 — 3. However, as in the pure rrih = 125 GeV case, a substantial 
enhancement of Rg^i'f'f) is most often associated with Rg+ A (ZZ) > Rg^i'fj) (contrary to the LHC observations). 
But this is not always the case. Among the rrih — itia scenarios we find 56 points in our parameter space for which 
Rg g l ' A (ZZ) < 1.3 and Rgg^ A (jj) > 1-3. Unfortunately for all those points the bb signal is predicted to be too strong, 
Rg+ A (bb) > 3.82. This situation is illustrated in Fig. [8j As seen from the upper panels in Fig. [8j for tan/3 = 1 there 
exist points (blue diamonds) such that r/, > 1 and R^J A {ll) > 1 ( or even > 1.5). However, the lower left panel of 
Fig. [8] shows that the Rg+ A (bb) values that correspond to those points are greater than 3.5, a result that is disfavored 
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TABLE VII: Table of maximum -Rg<^ A (77) values for the Type I 2HDM with m h = m,A = 125 GeV and associated R values 
for other initial and/or final states. The input parameters that give the maximal -Rg<J~ A (77) value are also tabulated. 
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TABLE VIII: Table of maximum #^(77) values for the Type II 2HDM with m h = m A = 125 GeV and associated R values 
for other initial and/or final states. The input parameters that give the maximal -Rg<J~ A (77) value are also tabulated. 
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FIG. 6: R Bg ^ A maximum values when = xua = 125 GeV as a function of tan/3 after imposing various constraints 
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see figure legend. Corresponding R Bg (ZZ) and R Bg (bb) are shown in the middle and lower panels 



by the LHC data. 

The case with wia ~ 125 GeV and mj? = 125 GeV is less attractive. For the Type I model, the constraints are 
such that once parameters are chosen so that H and A have the same mass of 125 GeV the maximum value achieved 
for -R^^ x (77) is rather modest reaching only 1.04 at small tan/3. For the Type II model, as seen in Fig. pi there are 
no parameter choices for which the H and A have a mass of 125 GeV while all other constraints are satisfied. 



C. rrih = 125 GeV and m,H = 125 GeV scenario 



Finally, we have the case where the h and H have a common mass of 125 GeV and we allow m#± and uia to vary 
freely. Following a similar search strategy, we find that some of the previously available tan/3 values are ruled out by 
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FIG. 7: As in Fig. [6] but for m„ = m A = 125 GeV. 
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the full set of constraints and that there is no gain in maximal Rg^ H ill) values, and often some loss, relative to the 
cases where only the h or only the H was required to have mass of 125 GeV. 



As discusssed earlier, the charged Higgs contribution to the 77 coupling loops is sometimes relevant. Therefore, in 
Fig. [9] we show separately the fermionic loop, W loop and H ± loop contributions normalized to the total amplitude 
for the most interesting cases of a Type I model with rrih = 125 GeV and with rrih — rriA = 125 GeV (left plots). One 
sees that the tan/3 values of 4 and 20 associated with Rg g (jj) ~ 1.3 are associated with large A H ± /A. Indeed, in these 
two cases, the relative charged Higgs contribution reaches nearly ~ 0.2 and is as large as the fermionic contribution, 
but of the opposite sign. In fact, although the dominant loop is the W loop, the loop may contribute as much as 
the dominant (top quark) fermionic loop. 

This should be contrasted with other cases, such as the Type II rrih — 125 GeV and rrih = rriA = 125 GeV cases 
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FIG. 8: Correlations between Rgj A {yi)/Rg g ( zz ) and Rgg (77) (upper panels) and RgJ (bb) (lower panels) for selected 
values of tan/?. 



illustrated in the right-hand plots of Fig. [9] One finds that the charged Higgs contributions are small when SUP 
constraints are imposed. In fact, the enhancement of Rg g {"fj) observed in Fig. fll prior to imposing SUP is caused just 
by the charged Higgs loop. When SUP constraints are imposed the charged Higgs amplitude is strongly reduced by 
the requirement that the quartic couplings not violate the perturbativity condition. Note that the SUP constraints 
can be violated even though all the mass parameters have been varied within what, a priori, appears to be a reasonable 
range, namely from a few GeV up to 1000 GeV. This is due to the fact that, for our input, the SUP conditions imply 
a strong constraint on m\ 2 that comes mainly from the requirement of keeping Ai small enough. 



IV. CONCLUSIONS 



We have analyzed the Type I and Type II two-Higgs-doublet extensions of the Standard Model with regard to 
consistency with a significant enhancement of the gluon-fusion-induced 77 signal observed at the LHC at ~ 125 GeV. 
All possible theoretical and experimental constraints have been imposed. We find that vacuum stability, unitarity and 
perturbativity play the key role in limiting the maximal possible enhancement which, in the most interesting scenarios, 
is generated by the charged Higgs loop contribution to the Higgs to two photon decay amplitude. Generically, we 
conclude that the Type II model allows a maximal enhancement of order of 2 — 3, whereas within the Type I model 
the maximal enhancement is limited to < 1.3. Current public results suggest an enhancement for gg — > h — > 77 that 
exceeds the maximum for the Type I model and would appear to favor the Type II model. 

However, we find that in the Type II model the parameters that give (-yy) > 1.3 are characterized by R^ g (ZZ) > 
Rg 0(77), a result that is inconsistent with experimental results in the gg — > h — > ZZ — > 4£ channel. Thus, if 
Rg g (-y^f) > 1.3 and R l g l g (ZZ) < 1.3 both persist experimentally, the Type II model cannot describe the data if only 
the h resides at 125 GeV. Similar statements apply to the case of the heavier H having a mass of 125 GeV. Next, 
we considered Type II models with approximately degenerate Higgs bosons at 125 GeV. We found that there exist 
theoretically consistent parameter choices for Type II models for which -R^ 1- " 4 (77) > 1.3 while R g l g ^ A {ZZ) < 1.3, but in 
these cases Rg^ A (bb) > 3.75, a value far above that observed. Thus, the Type II 2HDMs cannot yield R g l f A {jj) > 1.3 
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FIG. 9: For the most interesting scenarios we show imaginary part of charged Higgs contributions to the 77 amplitude normalized 
to the imaginary part of the sum of all (fermions, W + W~ , H + H~) contributions as a function of tan/3 after imposing all 
constraints. The parameters adopted correspond to maximal ^^(77) (or an appropriate sum for degenerate cases). 



without conflicting with other observables. In short, the Type II model is unable to give a significantly enhanced 
gg — > h — > 77 signal while maintaining consistency with other channels. 

In the case of the Type I model, the maximal R^ g {n) is of order of 1.3, as found if tan/3 = 4 or 20. In these 
cases, Rg g (ZZ) and R^ g {bb) are of order 1 as fairly consistent with current data. For these scenarios, the charged 
Higgs is light, m H ± — 90 GeV. Despite this small mass, there is no conflict with LHC data due to the fact that 
BR(t — > H + b) ~ 1/tan 2 (3 is small enough to be below current limits. Thus, Type I models could provide a consistent 
picture if the LHC results converge to only a modest enhancement for #^(77) < 1.3. But, if Rg g (j^f) is definitively 
measured to have a value much above 1.3 while the ZZ and bb channels show little enhancement then there is no 
consistent 2HDM description. One must go beyond the 2HDM to include new physics such as supersymmetry. 
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